Phase transitions of the binary production 2A — ► 3A, AA — > model 
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Phase transitions of the 2A —* 3A, AA — > reaction-diffusion model is explored by dynamical, 
iV-cluster approximations and by simulations. The model exhibits site occupation restriction and 
explicit diffusion of isolated particles. While the site mean-field approximation shows a single tran- 
sition at zero branching rate introduced in PRE 67, 056f 14 (2003), N > 2 cluster approximations 
predict the appearance of another transition line for weak diffusion (D) as well. The latter phase 
transition is continuous, occurs at finite branching rate and exhibits different scaling behavior. I 
show that the universal behavior of these transitions is in agreement with that of the PCPD model 
both on the mean-field level and in one dimension. Therefore this model exhibiting annihilation by 
quadruplets does not fit in the recently suggested classification of universality classes of absorbing 
state transitions in one dimension (PRL 90, 125701 (2003)). For high diffusion rates the effective 
2A —> 3A — > AA —* reaction becomes irrelevant and the model exhibits a mean-field transition 
only. The two regions are separated by a non-trivial critical endpoint at D* . 



Phase transitions in nonequilibrium systems, which do 
not possess hermitian Hamiltonian may appear in mod- 
els of population, epidemics, catalysis, cooperative trans- 
port [1], enzyme biology [2] and markets for example 
[3]. Reaction-diffusion systems are of primary interest 
since other nonequilibrium models often can be mapped 
onto them [4]. The classification of universality classes 
in reaction-diffusion systems [5,6] has recently got some 
impetus. In these systems particle creation, annihila- 
tion and diffusion processes compete and by tuning the 
control parameters phase transition may occur from an 
active steady state to an inactive, absorbing state of 
zero density. The fluctuations in the absorbing state 
are so small, that systems cannot escape from it, hence 
such phase transitions may emerge in one dimension al- 
ready. Several systems with binary, triplet or quadru- 
plet, particle reactions have been investigated numeri- 
cally and unclassified type of critical phase transitions 
were found [7-25]. Solid field theoretical treatment ex- 
ists for bosonic, binary production systems only [26], but 
this is not applicable for the active and critical states of 
site restricted models, since it cannot describe a steady 
state with finite density. 

The mean-field solution of general, 



iA A (n + k)A, 



iA — > (m - l)A, 



(1) 



models (with n > 1, m > 1, k > 0, I > and m — I > 0) 
resulted in a series of universality classes depending on n 
and m [24] . In particular for the n = m symmetrical case 
the density of particles above the critical point (a c > 0) 
scales as 



with j3 



MF 



p^\a-a c f, (2) 

1, while at the critical point it decays as 

P oc t~ a , (3) 

with a MF = i3 MF jvf F = l/n [23,24] (here "MF" de- 
notes mean- field value). On the other hand for the n < m 



asymmetric case continuous phase transitions at zero 
branching rate a c — occur with 



(3 MF = l/( m -n), a 



MF 



l/(m- 1) 



(4) 



For n > m the mean-field solution provides first order 
transition. 

The upper critical dimension for such systems is de- 
bated [6,21,23,24] but should be quite low (4 = 1- 2) 
allowing a few anomalous critical transitions only. For 
example d c < 1 was confirmed by simulations in case of 
the asymmetric, binary production 2A — > AA, AA — > 2A 
model [25]. It was also pointed out there that N > 1 
cluster mean-field approximation, that takes into account 
the diffusion of particles would provide a more adequate 
description of such models. Earlier studies have shown 
[28-30] that there exist models with first order transi- 
tions in the site mean-field approximation that changes 
to continuous one in higher level of cluster approxima- 
tions. Dependence on the diffusion was found to be im- 
portant in binary production models [8,10] and it turned 
out that at least N > 2 level of approximation is needed 
for an adequate description [17,20]. 

In this paper I investigate the 2A — > 3A, AA — > model 
and show that the diffusion plays an important role: it 
introduces a different critical point besides the one at 
a = branching rate. I show by simulations that this 
transition is not mean-field type in one dimension but 
belongs to the class of the: 2A — > 3A, 2A — > so called 
PCPD model. My model is defined and parametrized 
following the notation of [8] by the rules 



AA9, %AA -> AAA 
AAAA -► 0000 
A® » %A 



with rate a = (1 - p)(l - D)/2 

with rate A = p(A — D) 

with rate D . (5) 



Here D denotes the diffusion probability and p is the 
other control parameter of the system. 

Dynamical cluster mean-field approximations have 
been introduced for nonequilibrium models by [31,32]. 
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The master equations for N = 1, 2, 3, 4, 5 block probabil- 
ities were set up 

9 -^f^ = f(P»({ Si })), (6) 

where site variables may take values: Si — 0, A. Taking 
into account spatial reflection symmetries of P/v({sj}) 
this involves 20 independent variables in case of N = 5. 
The master equation (6) was solved numerically using the 
Runge-Kutta algorithm for N = 2,3,4,5 by several D 
and p values. The particle (p(p, D)) and pair {p2(p, D)) 
densities were determined by Pn({si})- For strong dif- 
fusion rates only a mean-field phase transition occurs at 
(7 = with (3=1/2 and a = 1/3 exponents belonging to 
the set set of classes (4) discovered in [24] . 

However for N > 1 and weak diffusion rates other 
phase transitions points emerge as well, with a c > 0. 
This means that for intermediate a and small D val- 
ues the absorbing state becomes stable as one can see 
in Fig.l. Simulations in one dimension confirm this (see 
later). 
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FIG. 1. Steady state density in N = 5 level approxima- 
tion for diffusion rates D = 0.5, 0.4, 0.35, 0.3, 0.2, 0.1, 0.05, 
0.01, 0.01 (top to bottom curves). Similar scenario appears 
for TV = 2, 3, 4. 



In the active phases in the neighborhood of the a c > 
transition points power law fitting of the form (2) to the 
mean-field data point resulted in (3 = 1 for all N > 1 
levels of approximations. On the other hand for the pair 
density in pair approximations one obtains (3=1 again, 
like in case of the PCPD model for low diffusions [8]. 
This anomaly disappears for N = 3, 4, 5 and the fitting 
results in (3 = 2 for pairs. 

At the a c > critical points the dynamical behav- 
ior is power-law type (3) with a = 1/2 for N = 3,4,5. 
Again the pair approximation gives the strange result: 
a = 1 (like in ref. [8]). The failure of the pair approxima- 
tion also appears in the inactive region, where it results 
in exponential density decay. In contrast with this the 



N = 3, 4, 5 approximations show power-laws here with 
a = 1 for particles and a = 2 for pairs. The above 
a and (3 exponents occurring by low diffusions at the 
critical points and inactive phases in the phase diagram 
away from the a = transition are different from those 
of the site mean-field values (4) . This can be explained 
by accepting that the dominant decay process for a > 0, 
D < D* is: 2A -> (via 2A -> 3A -> AA -> 0) instead 
of the A A — + - that is the only mode of decay at a = 0. 
Altogether one can find very similar cluster mean-field 
behavior as in case of the PCPD model [8,20]. 

One can also observe that by increasing D from zero 
the PCPD like transitions disappear at some D* value, 
when the p(oo) steady state curve touches the p = axis. 
For D > D* there is no absorbing state in the system and 
a critical endpoint appears with (3 = 2 (parabolic) sin- 
gularity at a*. For N = 5 the endpoint is located at 
D* = 0.301(1), =0.53(1). 

To test these analytical findings I have performed sim- 
ulations in one dimension. These were carried out on 
L = 1 — 5 x 10 5 sized systems with periodic boundary 
conditions. The initial states were half filled lattices with 
randomly distributed As and the density of particles is 
followed up to 5 x 10 8 Monte Carlo steps (MCS). One 
MCS consist of the following processes. A particle, a di- 
rection and a number x G (0,1) are selected randomly; 
if x < D a site exchange is attempted with one of the 
randomly selected empty nearest neighbors (nn); else if 
D < x < (D + A) four neighboring particles are removed; 
else one particle is created at an empty site in the ran- 
domly selected direction following a pair of As. In each 
MCS the time is updated by 1/n, where n is the number 
of particles. 




FIG. 2. cteff in the one dimensional 3A — > iA, 4A — » 
model at D = 0.5 Different curves correspond to p = 0.1583, 
0.1584, 0.1585, 0.15852, 0.15853, 0.1586, 0.1587, 0.159 (top to 
bottom) 

First I followed the density of particles for a small a 
(at p = 0.95) at diffusion rates: D = 0.5 and D = 0.2. 
In both cases a power-law decay with a = 0.5 exponent 
could be observed, hence an inactive phase with decay 
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of the AA — > process — valid in one dimension [34] - 
was identified. 

The critical points were determined by calculating the 
local slopes defined as 



«//(*) 



\n[p{t)/p{t/m)} 
ln(m) 



(J) 



(where I used m = 2) for D = 0.2, 0.5, 0.747. As Fig.2 
shows the local slopes curve for D = 0.5, p = 0.15850(2) 
extrapolates to a = 0.21(1). This value agrees with 
that of the PCPD model [20,21]. Other curves exhibit 
curvature for long times, i.e. for p < 0.1585 they veer 
up (active phase), while for p > 0.1585 they veer down 
(absorbing phase). The local slopes figure shows simi- 
lar strong correction to scaling as in case of the PCPD 
model, i.e. some curves that seem to be supercritical veer 
down after t >~ 10 6 MCS. Similar results are obtained 
by other a c > transitions. For D = 0.2, when the 
critical point is at p — 0.0892(1) the local slopes for the 
density decay predicts a = 0.21(2). Repeating the sim- 
ulations at D = 0.9 no absorbing phase has been found 
(up to p < 0.9999), the steady state density disappears 
monotonously as a —> 0. At a = the density decays 
with a = 1/3 valid for the 4 A — > process in one dimen- 
sion [35]. 

The steady state density in the active phase near the 
critical phase transition point is expected to scale as 



p(oo) 



oc \p - p t 



1/3 



Using the local slopes method one 



can get a precise estimate for j3 and see the corrections 
to scaling 



PeffiPi) = 



ln.p(oo,pi) - lnp(oo,p t _i) 
ln(p 4 ) - ln(pi_i) 



(8) 



The steady state density was determined by running the 
simulations in the active phase: e — p c — pi > 0, by aver- 
aging over ~ 100 samples in a time window following the 
level-off is achieved. As one can see on Fig. 4 the effective 
exponent tends to = 0.40(2) as e -> both for D = 0.5 
and D — 0.2 diffusions. These values are in agreement 
with that of the one dimensional PCPD model [20,21]. 
Again assuming logarithmic corrections as in [20] of the 
form 



p(oo,e) = [e/(a + bln(e))f 



(9) 



one can obtain j? = 0.1585(1) and/3 = 0.38(1) for D = 0.5 
and p = 0.0892(1) and fi = 0.41(3) for D = 0.2, which 
agrees with the previous values within numerical accu- 
racy. [11,24]. Altogether one can not see relevant log- 
arithmic corrections for the diffusion rates investigated 
here. 

In case of D = 0.9, <r c = one can see {3 = 0.50(1) 
in agreement with the N = 3, 4, 5 cluster mean-field ap- 
proximation results. A quadratic fitting of the form 



results in: a = 0.195, b = 0.158, x = 0.214, (3 = 0.51(1). 
This suggests that the effective 2A — > process is weaker 
now than the 2 A — > 3 A, leaving the transition at <j c = 0. 
The phase diagram for different levels of approximations 
as well as MC data are shown on Fig. 3. As one can see 
approximations tend towards the simulated points by in- 
creasing N. 
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FIG. 3. Phase diagram. Stars correspond to N — 2, boxes 
to N = 3, bullets to N — 4 and triangles to N — 5 cluster 
mean-field approximations. Diamonds denote Id simulation 
data. The lines serve to guide the eye. At the a = line a 
mean-field transition occurs. 



Similar reentrant phase diagram has been observed in 
case of the unary production, triplet annihilation model 
(A -> 2A, 3A -> 0) [36] and in a variant of the NEKIM 
model [37]. In all cases the diffusion competes with par- 
ticle reaction processes, and the bare parameters should 
somehow form renormalizcd reaction rates which govern 
the evolution over long times and distances, the details 
have not been worked out. 




Peft = P - 



be 2 



(10) 



FIG. 4. /3 e ff as the function of e in the one dimensional 
2A -» 3^4, AA -> model. The bullets correspond to D = 0.5, 
the boxes to D = 0.2, the diamonds to D = 0.9 diffusion rate. 
The solid line shows a quadratic fitting of the form (10). 
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Finite size scaling investigations at D = 0.5 and 
p c = 0.1585 were performed for system sizes: Li = 
32, 64, 128. ..4096. The quasi-steady state density (aver- 
aged over surviving samples) is expected to scale accord- 
ing to 



p a (oo,p c ,L)<xL-^, 



(11) 



while the characteristic lifetime for half of the samples 
to reach the absorbing state scales with the dynamical 
exponent Z as 



T(p c ,L) oc L z . 
These quantities were analyzed by the local slopes: 
lnr(ij) - Inr(ij_i) 



Zeff(L) 
(3/v ± (L) 



hi Li — lnLj_i 
lnp s (L t ) - lnp s (£ t _i) 
In Li — In Lj_i 



(12) 

(13) 
(14) 



Linear extrapolation to L — > oo results in Z = 1.80(15) 
and /3/v± = 0.40(3). These values corroborate that the 
transition is of PCPD type 

In conclusion the N cluster mean-field study of the bi- 
nary production 2A — > 3A, AA — > model has shown the 
appearance of another critical transition with non-zero 
production rate for low diffusions. While the pair approx- 
imation results in somewhat odd results - like in case of 
other binary production systems - the N — 3, 4, 5 levels 
coherently exhibit PCPD-likc mean-field critical behavior 
for these phase transition points and within the absorb- 
ing phase. This transition line disappears at a critical 
endpoint for D > D* characterized by (3 — 2 order pa- 
rameter singularity and for high diffusion rates the <j c = 
critical point remains only in the system, predicted by the 
site mean-field approximation. The utmost importance 
of diffusion dependence and the corresponding N > 2 
cluster mean-field approximations is demonstrated in this 
study. 

Extensive simulations in one dimension have confirmed 
the existence of the nontrivial transition for low diffu- 
sions. By these transitions points the critical behavior 
agrees with that of the latest results obtained for the 
PCPD model. Therefore this model does not fit in the 
table of universality classes suggested for such models in 
one dimension [21]. The reason behind this discrepancy 
might be that in [21] low diffusions have not been inves- 
tigated or the lack of complete site exclusion in their 
model. Site exclusion has been shown to be relevant 
in multi-species reaction-diffusion systems and in binary 
production systems [38] . 

An interesting, open problem is the exploration of the 
phase structure of this system in higher dimensions. The 
agreement of one-dimensional results with those of the 



cluster mean-field shows that similar rich phase struc- 
ture may emerge in higher dimensions too. That would 
mean that the effective 2A — > reaction is generated 
via: 2A — > 3^4 — > AA — > again. These results raise the 
possibility that such mechanism also emerges by unary 
production systems (example by: A — ► 2A, AA — > 0) and 
one should find a DP transition instead of the mean-field 
one suggested by perturbative renormalization study [35] 
of such models. This would affect the classification of 
fundamental universality classes of RD systems and may 
point out the weak points of the perturbative renormal- 
ization. An other important point to be investigated is 
the scaling behavior at the critical endpoint. 
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